Geometric quantization of time-dependent completely integrable 
Hamiltonian systems 
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A time-dependent completely integrable Hamiltonian system is quantized with respect to 
Q\ | time-dependent action-angle variables near an instantly compact regular invariant manifold. 

Its Hamiltonian depends only on action variables, and has a time-independent countable 
■ energy spectrum. 
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A time-dependent Hamiltonian system of m degrees of freedom is called a completely 
integrable system (henceforth CIS) if it admits m independent first integrals in involution. 
Choosing appropriate dynamic variables, one may hope to quantize a time-dependent CIS 
Q_i| so that its quantum Hamiltonian and first integral operators possess time-independent 
spectra. 1 Time-dependent action-angle variables introduced below are of this type. Written 
relative to these variables, a Hamiltonian of a time-dependent CIS is a function only of the 
action coordinates. It follows that, if time- dependent action-angle coordinates hold fixed, a 
time-dependent CIS can be quantized just as an autonomous one, and its energy spectrum 
. !_h \ is time-independent. 

In order to introduce time-dependent action-angle variables, we use the fact that a 
time-dependent CIS of m degrees of freedom can be extended to an autonomous one of 
m + 1 degrees of freedom where the time is regarded as a dynamic variable. 2-4 By virtue 
of the classical Arnold-Liouville theorem, 5 ' 6 an autonomous CIS admits the action-angle 
coordinates around a regular connected compact invariant manifold. The problem is that 
invariant manifolds of a time-dependent CIS are not compact because of the time axis. 
Therefore, we first generalize the above mentioned theorem to noncompact invariant mani- 
folds. Then we show that, if a regular connected invariant manifold N of a time- dependent 
CIS is compact at each instant, it admits an open neighbourhood in the ambient momentum 
phase space which is isomorphic to the product 

^ = 1 x T m x ^ (1) 



Electronic mail: fiorani@mat.unimi.it 
Electronic mail: giachetta@campus.unicam.it 

Electronic mail: sard@grav.phys.msu.su; URL: http://webcenter.ru/^sardan/ 



1 



of the time axis M, an m-dimensional torus T m and an open domain V C M. m . This product 
is equipped with the coordinates 

li), i = l,...,m, (2) 

where t is the Cartesian coordinate on M and tp 1 are cyclic coordinates on T m . Written with 
respect to these coordinates, a Hamiltonian and first integrals of a time- dependent CIS are 
functions only of action coordinates l{. The corresponding Hamilton equation on W reads 

k = 0, jj = &H{Ij). 

A glance at this equation shows shows that, given action-angle coordinates (Q), a time- 
dependent CIS can be seen as an autonomous CIS on the symplectic annulus 

p = V x T m , (3) 

equipped with the action-angle coordinates Jj) and provided with the symplectic form 

VI P = dli A d<f>*. (4) 

Therefore, we can quantize a time-dependent CIS with respect to action-angle variables 
similarly to that of an autonomous CIS. Of course, the choice of time-dependent action- 
angle coordinates by no means is unique. They are changed by canonical transformations. 
Therefore, we employ the geometric quantization technique 7-9 which remains equivalent 
under such kind transformations. At the same time, geometric quantization essentially 
depends on the choice of polarization. 10,11 

Geometric quantization of an autonomous CIS has been studied with respect to polar- 
ization spanned by Hamiltonian vector fields of first integrals. 12 In fact, the Simms quanti- 
zation of the harmonic oscillator 9 is also of this type. The problem is that the associated 
quantum algebra includes functions which are not defined on the whole momentum phase 
space, and elements of the carrier space fail to be smooth sections of the quantum bundle. 
Indeed, written with respect to the action-angle variables, this quantum algebra consists of 
functions which are affine in angle coordinates. 

We choose a different polarization spanned by almost-Hamiltonian vector fields d k of 
angle variables. The associated quantum algebra A consists of smooth functions which are 
affine in action variables. Note that this quantization of the symplectic annulus P (|3|) is 
equivalent to geometric quantization of the cotangent bundle T*T m of the torus T m with 
respect to the familiar vertical polarization. As is well-known, the vertical polarization of 
a cotangent bundle leads to its Schrodinger quantization. We show that A possesses a 
set of nonequivalent representations in the separable pre-Hilbert space C°°(T m ) of smooth 
complex functions on T m . In particular, the action operators read 

I k = -id k - A fc , (5) 
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where are real numbers which specify different representations of A. By virtue of the 
multidimensional Fourier theorem, an orthonormal basis for C°°(T m ) consists of functions 

V>( nr ) = exp[i(n r (j) r )], (n r ) = (m, . . . , n m ) G IT 1 . (6) 

With respect to this basis, the action operators @ are countable diagonal matrices 

4^V) = ( nfc ~ ^)VV)- ( 7 ) 

Given the representation (|5|), any polynomial Hamiltonian 7i{Ik) of a CIS is uniquely quan- 
tized as a Hermitian element H(h) = H.{Ik) of the enveloping algebra A of A. It has the 
countable time-independent spectrum 

^(4)VV) = ^(nr)VV)> ^(nr) = - A fe ), n fc G (ra r ). (8) 

Similarly, polynomial first integrals are quantized. Since Ik are diagonal, one can also 
quantize Hamiltonians Ti.(Ij) and first integrals F(Ij) which are analytic functions on R m . 

Note that, because geometric quantization is equivalent under canonical transforma- 
tions, quantization of a time-dependent CIS with respect to action-angle variables induces 
its quantization with respect to initial variables near an invariant manifold in the ambient 
momentum phase space. However, its Hamiltonian need not be represented in terms of first 
integrals and need not belong to the quantum algebra A because it fails to be a scalar under 
time-dependent transformations. Moreover, this induced quantization can not be in general 
extended to the whole momentum phase space because of the topological obstructions to 
the existence of global action-angle coordinates. 13,14 

For instance, one usually mentions a harmonic oscillator as the simplest CIS whose 
quantization in the action- angle variables looks notoriously difficult. 15 However, a harmonic 
oscillator written relative to action-angle coordinates ((f), 1) is located in the momentum 
phase space R 2 \ {0}, but it is not the standard oscillator on R 2 . Namely, there is a 
monomorphism, but not an isomorphism of the Poisson algebra of smooth complex functions 
on R 2 to that on R 2 \{0}. In particular, the angle polarization on R 2 \{0} is not extended to 
R 2 . As a consequence, the quantum algebra associated to this polarization is not extended 
to R 2 , and so is its carrier space C°°(T m ). 

In conclusion, let us remark that, since Hamiltonians depend only on action variables and 
possess time-independent countable spectra, quantum CISs look especially promising for 
holonomic quantum computation, based on driving degenerate eigenstates of a Hamiltonian 
over a parameter space. 16-18 We will construct the corresponding quantum control operator. 

II. CLASSICAL COMPLETELY INTEGRABLE SYSTEMS 

Recall that the configuration space of time-dependent mechanics is a fiber bundle Q — > R 
over the time axis R equipped with the bundle coordinates (t,q k ), k = l,...,m. The 
corresponding momentum phase space is the vertical cotangent bundle V*Q of Q — > R 
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endowed with holonomic coordinates (t, q k ,p k ). 19 ' 20 The cotangent bundle T*Q, coordinated 
by (q x ,p\) = (t,q k ,p ,p k ), plays a role of the homogeneous momentum phase space. It is 
provided with the canonical Liouville form H = p\dq x , the symplectic form Q = oE, and 
the corresponding Poisson bracket 

{/, f'h = d x fd x f - d x fd x f, f, f G C°°(T*Q). (9) 

There is the one- dimensional trivial affine bundle 

C : T*Q -> V*Q. (10) 

Given its global section h, one can equip T*Q with the global bundle coordinate r = p — h. 

The fiber bundle ( |10"D provides the vertical cotangent bundle V*Q with the canonical 
Poisson structure {, }y such that 

C{f, f'}v = {Cf, Cfh, V/, /' G C°°(V*Q), (11) 
{fj'}v = d k fd k f'-d k fd k f. (12) 

Its characteristic symplectic foliation coincides with the fibration V*Q — > R. However, the 
Poisson structure (p!2|) fails to set any dynamic equation on the momentum phase space 
V*Q because Hamiltonian vector fields 

•&j = d k fd k - d k fd k , # f \df = {f, f} v , f, f g c°°(v*g), 

of functions on V*Q are vertical. 

A Hamiltonian of time- dependent mechanics is defined as a global section 

h : V*Q - T*Q, Po oh = -Hit, q j , Pj ), 

of the affine bundle ( (0)- 19,20 It yields the pull-back Hamiltonian form 

H = h*E = p k dq k - Hdt (13) 

on V*Q. Then there exists a unique vector field •jh on V*Q such that 

j H \dt = l, j H \dH = 0, 

lH = d t + d k Hd k -d k Hd k . (14) 
Its trajectories obey the Hamilton equation 

q k = a k n, p k = -d k n. (is) 



A first integral of the Hamilton equation ( p!q) is a smooth real function F on V*Q whose 
Lie derivative 



L lH F = lH \dF = d t F + {H,F} 



v 
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along the vector field (0) vanishes, i.e., F is constant on trajectories of 7#. A time- 
dependent Hamiltonian system (V*Q, H) is said to be completely integrable if the Hamilton 
equation flI5| ) admits m first integrals F k which are in involution with respect to the Poisson 



bracket {, }y (|I2 ) and whose differentials dF k are linearly independent almost everywhere, 
i.e., the set of points where this condition fails is nowhere dense in V*Q. One can associate 
to this system an autonomous CIS on T*Q as follows. 



Let us consider the pull-back (*H of the Hamiltonian form H (|13|) onto the cotangent 
bundle T*Q. It is readily observed that 

H* = d t \(E-Ch*E)=p + n (16) 

is a function on T*Q. Let us regard 7i* as a Hamiltonian of an autonomous Hamiltonian 
system on the symplectic manifold (T*Q,Q). Its Hamiltonian vector field 

7 t = d t - d t Hd° + d k Hd k - d k Hd k (17) 



is projected onto the vector field 7^ (|14]) on V*Q so that 

C*(W) = {H\ Cfh, V/ e C°°{V*Q). 
An immediate consequence of this relation is the following. 

Proposition 1: (i) Given a time-dependent CIS (V*Q, H; F k ) on V*Q, the Hamiltonian 
system (T*Q; TC*, (*F k ) on T*Q is completely integrable. (ii) Let N be a connected regular 
invariant manifold of (V*Q, H; F k ). Then h(N) C T*Q is a connected regular invariant 
manifold of the autonomous CIS (T*Q;H*,(*Fk)- 

Hereafter, the vector field 7^ (|14|) is assumed to be complete. In this case, the Hamilton 
equation (|TJ|) admits a unique global solution through each point of the momentum phase 
space V*Q, and trajectories of 7# define a trivial bundle V*Q — > V t *Q over any fiber V t *Q 
of V*Q — ► R. Without loss of generality, one can choose the fiber i : VqQ — > V*Q at t = 0. 
Since iV is an invariant manifold, the fibration 

£ : V*Q -> V *g (18) 

also yields the fibration of iV onto iV = A^fl V$Q such that = R x AT is a trivial bundle. 

III. TIME-DEPENDENT ACTION-ANGLE COORDINATES 



Let us introduce the action-angle coordinates around an invariant manifold N of a time- 
dependent CIS on V*Q by use of the action-angle coordinates around the invariant manifold 
h(N) of the autonomous CIS on T*Q in Proposition |l| Since N and, consequently, h(N) 
are noncompact, we first prove the following. 

Proposition 2: Let M be a connected invariant manifold of an autonomous CIS {Fx}, 
A = 1, . . . , n, on a symplectic manifold (Z, flz), and let the Hamiltonian vector fields of the 
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first integrals F a on M be complete. Let U be a neighbourhood of M such that {F\} have 
no critical points in U and the submersion xF\ : U — > M™ is a trivial bundle of Lagrangian 
invariant manifolds over a domain V' C M n . Then U is isomorphic to the symplectic annulus 

W' = R n ~ m xT m xV (19) 

provided with the generalized action-angle coordinates 

(x\...,X n - m ,r- m+ \...,<p n , h,...,In) (20) 

such that the symplectic form on W reads 

Vt z = dl a A dx a + dl k A d(p h , 
and the first integrals F\ depend only on the action coordinates (/«). 

Proof: In accordance with the well-known theorem, 5,21 the invariant manifold M is diffeo- 
morphic to the product IR™"™ x T m , provided with coordinates (y x ) = (s a , if 1 ) where ip l are 
linear functions of parameters s x along the integral curves of Hamiltonian vector fields of 
first integrals F\ on M. Let (J\) be coordinates on V which are values of first integrals 
F x . Since W — ► V' is a trivial bundle, (y x , J\) are bundle coordinates on the annulus W 
(PH)- It should be emphasized that, since group parameters are given up to a shift, the 
coordinates y x on W are determined up to a shift by functions of coordinates J\. Written 
relative to these coordinates, the symplectic form Qz on W reads 

Vt z = n af3 dJ a A dJ p + Q%dJ a A dyP. (21) 

By the definition of coordinates (y x ) , the Hamiltonian vector fields i)\ of first integrals take 
the coordinate form $ A = i?"(J M )9 Q where 

#a = d a + i a (J x )d i , ^ = ^(Jx)d k , (22) 

and they obey the relations 

$ x \tt z = -dJx, fi|0? = ^. (23) 

It follows that VL'p is a nondegenerate matrix and = (fi^ 1 )^, i.e., the matrix functions 
f2jg depend only on coordinates J\. In Appendix A, we obtain the desired coordinates 

X a = s\ <j) a (s b ^\Jx), I a = Ja, Wk)- (24) 

Note that, if M is a compact invariant manifold, the conditions of Proposition |2| always 
hold. 6 
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Of course, the generalized action-angle coordinates fl2"0| ) by no means are unique. For 
instance, let JF a , a = 1, . . . , n — m be an arbitrary smooth function on R m . Let us consider 
the canonical coordinate transformation 

x la = x\ 0' J = J + x a ^ a (/,), i£ = 1.-^(1,-), 4 = / fc . (25) 

Then (x /a , J„, JjQ also generalized action-angle coordinates on the symplectic annulus 
which differs from W ( jl9|) in another trivialization. 

Now, we apply Proposition |2] to the CISs in Proposition |l[ 

Proposition 3: Let N be a connected regular invariant manifold of a time-dependent CIS 
(V*Q, H; Fk), and let the image N of its projection £ ([1S|) be compact. Then the invari- 
ant manifold fo(iV) of the autonomous CIS (T*Q;TC* , C*Jfc) nas an open neighbourhood U 
obeying the condition of Proposition [| 

The proof is in Appendix B. In accordance with Proposition the open neighbourhood 
U of the invariant manifold h(N) in Proposition ^ is isomorphic to the symplectic annulus 

W' = R x T m x V', (26) 

provided with the generalized action-angle coordinates (t, 1 , . . . , m , Jo, ■ ■ ■ , J m ) such that 
the symplectic form on W reads 

Q! = dl A alt + difc A rf0 fc . 

A glance at the Hamiltonian vector field $0 = 7t ([17|) and the relation (|23|) - (0), shows 
that Iq = Jo = 7~C* and the corresponding generalized angle coordinate is x° = t, while the 
first integrals = (*Fk depend only on the action coordinates Jj. 

Since the action coordinates Jj are independent of the coordinate J , the symplectic 
annulus W fl2l)| ) inherits the fibration (|TDp which reads 

C:W'3 (t,0\ J , J*) -> e W = R x T m x V. 

By the relation similar to (ITTJ), the product W is provided with the Poisson structure 

{/, /V = - djd'f, /, /' e C°°(W). (27) 

Therefore, one can regard W with coordinates (t, 0*, Jj) as the momentum phase space of 
the time-dependent CIS in question around its invariant manifold N. 

It is readily observed that the Hamiltonian vector field 7^ of the autonomous Hamilto- 
nian TC* = Io is 7t = d t , and so is its projection 7^ (|TJ|) on H 7 . Consequently, the Hamilton 
equation (p~5f ) with respect to the action- angle coordinates take the form Jj = 0, <$> % — 0. 
Hence, (i, <j) 1 , Jj) are the initial date coordinates. One can introduce such coordinates as 
follows. Given the fibration £ (|T8|), let us provide N x V C Vq*Q in Proposition [| with 
action-angle coordinates ((f) , Jj) for the CIS {ig F^} on the symplectic leaf V *Q. Then, it is 
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readily observed that (t, 4> , U) are time-dependent action-angle coordinates on W (P such 
that the Hamiltonian T~C{Ij) of a time-dependent CIS relative to these coordinates vanishes, 
i.e., 7i* = Iq. Using the canonical transformations (0), one can consider time-dependent 
action-angle coordinates besides the initial date ones. Given a smooth function 7i on M. m , 
let us further provide W with the action-angle coordinates 

p = +t&H<J j ), J = 7 -W(7 i ), k = h 

such that TC(Ii) is the Hamiltonian. 

Thus, action-angle coordinates for a time-dependent CIS provide a particular solution of 
the problem of a representation of its Hamiltonian in terms of first integrals. 22 ' 23 However, 
this representation need not hold with respect to the initial bundle coordinates on V*Q 
because a Hamiltonian fails to be a scalar under time- dependent transformations. 

IV. QUANTUM COMPLETELY INTEGRABLE SYSTEMS 

In order to quantize a time-dependent CIS on the Poisson manifold (W, {,}w), one 
may follow the general procedure of instantwise geometric quantization of time- dependent 
Hamiltonian systems in Ref. 24. As was mentioned above, it however can be quantized 
as an autonomous CIS on the symplectic annulus (P,Q P ) with respect to fixed time- 
dependent action-angle coordinates. 

In accordance with the standard geometric quantization procedure, 7 ' 8 since the sym- 
plectic form Qp (|) is exact, the prequantum bundle is defined as a trivial complex line 
bundle C over P. Since the action-angle coordinates are canonical for the symplectic form 
(£D, the prequantum bundle C need no metaplectic correction. Let its trivialization 

C = P x C (28) 

hold fixed. Any other trivialization leads to equivalent quantization of P. Given the 
associated bundle coordinates ((p k ,Ik,c), c G C, on C ([28|), one can treat its sections as 
smooth complex functions on P. 

The Konstant-Souriau prequantization formula associates to each smooth real function 
/ G C°°(P) on P the first order differential operator 

f = -iV* f + f (29) 

on sections of C, where = d k fdk — dkfd h is the Hamiltonian vector field of / and V is 
the covariant differential with respect to a suitable ?7(l)-principal connection on C. This 
connection preserves the Hermitian metric g(c, d) = cd on C, and its curvature form obeys 
the prequantization condition R = iQp. It reads 

A = A + icl k d(j) k <g> d c , (30) 
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where A is a flat [/(l)-principal connection on C — > P. The equivalence classes of flat 
principal connections on C are indexed by the set M m /Z m of homomorphisms of the de 
Rham cohomology group 

H\P) = H\T m ) = R m 

of P to the cycle group U(l). 9 We choose their representatives of the form 

A)[(Afc)] = dl k ®d k + d<f) k <g> + iA fc c<9 c ), A fc g [0, 1). 

Then the connection (j30|) up to gauge transformations reads 

A[(X k )\ = dl k <g> d k + d<\> k ® + i{I k + A fc )c«9 c ). (31) 

For the sake of simplicity, we will assume that the numbers A^ in the expression belong 
to E, but will bear in mind that connections A[(Afe)] and A[(A^.)] with \ k — X' k G Z are 
gauge conjugated. Given a connection Q5ID, the prequantization operators read 

/ = -^ / + (/-(4 + A fc )9 fc /)- (32) 

Let us choose the above mentioned angle polarization Vir which is the vertical tangent 
bundle of the fibration 7r : P —>■ T m , and is spanned by the vectors d k . It is readily observed 
that the corresponding quantum algebra A C C°°(P) consists of affrne functions 

f = a k (^)I k + b(^) (33) 



of action coordinates I k - The carrier space of its representation by operators ( |32"D is defined 
as the space E of sections p of the prequantum bundle C of compact support which obey 
the condition V^p = for any Hamiltonian vector field $ subordinate to the distribution 
Vtt. This condition reads 

d k fd k p = 0, V/ G C°°(T m ). 

It follows that elements of E are independent of action variables and, consequently, fail to 
be of compact support, unless p = 0. This well-known problem of Schrodinger geometric 
quantization is solved as follows. 24,25 

Fix a slice %t '■ T m — > T m x V. Let Ct = i^C be the pull-back of the prequantum bundle 
C ( p8|) over the torus T m . It is a trivial complex line bundle Ct = T m x C provided with 
the pull-back Hermitian metric g(c, d) = cc'. Its sections are smooth complex functions on 
T m . Let 

A = i* T A = d(f) k <g> (dk + i(I k + A fc )c<9 c ) 

be the pull-back of the connection A ([$]]) onto Ct- Let Dbea metalinear bundle of complex 
half-forms on the torus T m . It admits the canonical lift of any vector field r on T m , and 
the corresponding Lie derivative of its sections reads 
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Let us consider the tensor product 

Y = C T ®V^T m . (34) 

Since the Hamiltonian vector fields 

= a k d k - (I r d k a r + d k b)d k 



of functions / (|33|) are projectable onto T m , one can associate to each element / of the 
quantum algebra A the first order differential operator 

/ = (-iV^ + /) <g> Id + Id ® = -m fc <9 fc - ^<9 fc a fc - a k \ k + 6 (35) 



on sections of Y. A direct computation shows that the operators (|35|) obey the Dirac 
condition 

[/,/'] = -i{fj'}- 



Sections px of the quantum bundle Y — > T m (p4|) constitute a pre-Hilbert space -Sr with 
respect to the nondegenerate Hermitian form 



(pt\p't) = Q^j J PtPt, pt,p' t e #t- 



Then it is readily observed that / (|35D are Hermitian operators in E?- In particular, the 
action operators take the form (|5|). 

Of course, the above quantization depends on the choice of a connection A[(Afc)] ( |3~1~D 
and a metalinear bundle T>. The latter need not be trivial. If T> is trivial, sections of the 
quantum bundle Y — > T m ( |34"D obey the transformation rule 

p T {<p k + 2tt) = p r (0 fc ) 

for all indices k. They are naturally complex smooth functions on T' m . In this case, Ej 1 is 
the above mentioned pre-Hilbert space C°°(T m ) of complex smooth functions on T m whose 
basis consists of functions (^). The action operators / (^) with respect to this basis are 
represented by countable diagonal matrices (|7]), while functions a(4>) are decomposed into 
the pull-back functions ip(n r ) which act on C°°(T m ) by multiplications 

1p(n r )1p(n' r ) = VVrMn'r) = VV+O' (36) 



If T> is a nontrivial metalinear bundle, sections of the quantum bundle Y — > T m 
obey the transformation rule 

Pt(0 j +27t) = - Pt ^) (37) 
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for some indices j. In this case, the orthonormal basis of the pre-Hilbert space Ef can be 
represented by double-valued complex functions 



Aniw) = exp[i(ni0 l + (rij + )^W)} (38) 
on T m . They are eigenvectors 

of the operators Ik and the functions a(0) act on the basis (|38|) by the above law 
(^). It follows that the representation of A determined by the connection A[(Afc)] (|T]) 



in the space of sections ([J7J) of a nontrivial quantum bundle Y fl3"4|) is equivalent to its 
representation determined by the connection A[(Aj, Xj — |)] in the space C°°(T m ) of smooth 
complex functions on T m . 

Therefore, one can restrict the study of representations of the quantum algebra A to its 
representations in C°°(T m ) associated to different connections (|3"T|). These representations 
are nonequivalent, unless Xk — A^. G Z for all indices k. 

Now, in order to quantize the Poisson manifold (W, {,}w), one can simply replace 
functions on T m with those on R x T m . 7 ' 2A Let us choose the angle polarization of W 
spanned by the vectors d k . The corresponding quantum algebra Aw C C°°(W) consists of 
affine functions 

/ = a k (t^)I k + b(t^) (39) 



of action coordinates Ik, represented by the operators (p5|) in the space C°°(1R x T m ) of 
smooth complex functions on R x T m . This space is provided with the structure of the 
pre-Hilbert C°°(R)-module with respect to the nondegenerate C°°(R)-bilinear form 



m>) = (to) J ^ ' ^ ^ e C °° (M x T ' 



Its basis consists of the pull-backs onto R x T m of the functions ip( nr ) (§)• 

Since the Poisson structure ([27]) defines no dynamics on the momentum phase space 
W ([!]), we should quantize the homogeneous momentum phase space W (^) in order 
to describe evolution of a quantum time- dependent CIS. Following the general scheme in 
Refs. 25,26, one can provide the relevant geometric quantization of the symplectic annulus 
(W',Q'). The corresponding quantum algebra Aw' C C°°(W') consists of affine functions 

f = a\t^)h + b{t^) 

of action coordinates I\. It suffices to consider its subalgebra consisting of the elements / 
and Jo + / for all / £ Aw d39|). They are represented by the operators / ( p5|) and Iq = —id t 
in the pre-Hilbert module C°°(R x T m ). If a Hamiltonian TC(Ij) of the time-dependent 
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CIS is a polynomial (or analytic) function in action variables, the Hamiltonian 7i* of the 
associated autonomous CIS is quantized as 

H* = -id t + H(Tj). 

Then we obtain the Schrddinger equation 

H*i) = -id t i) + H(-id k - A fc )^ = 0, ^eC°°(lx T m ). 

Its solutions are the series 

V> = B M exp[-iE {nr) t]^j {nr) , B {nr) e C, 

(n r ) 

where Ef n \ are the eigenvalues (||) of the Hamiltonian 7i. 

In conclusion, bearing in mind applications to holonomic quantum computation, let 
us choose action-angle coordinates such that a Hamiltonian 7i of a CIS is independent of 
action variables I a (a,b,c = 1, . . . , /). Then its eigenvalues are countably degenerate. Let 
us consider the perturbed Hamiltonian 



3)1 



where the perturbation term A depends on the action-angle coordinates with the above 
mentioned indices a, b, c, . . . and on some time-dependent parameters s M (t) by the law 



A = A a Js^ ( j ) b )d t s^ a . 



(40) 



The Hamiltonian 7i' characterizes a CIS with time-dependent parameters. 20,26,27 Being 
affine in action variables, the perturbation term A ( f4"0"| ) is represented by the operator 

A = -{ik a p d a + % -d a A% + \ a K%)d t sK 

Since the operators A and 7i mutually commute, the total quantum evolution operator falls 
into the product 





t 




t 




t 


Texp 


-% J n'dt' 


= Texp 


-i J ndt' 


o Texp 


-i J Adt' 


















The first factor in this product is the dynamic evolution operator of the quantum CIS. The 
second one acts in the eigenspaces of the dynamic Hamiltonian n and reads 



Texp 



{-A a p (<P b , s»(t'))d a - -dA^, s"(f )) + i\ a A%(<f> b , s»{t'))}d t s?dt' 



Texp 



a([0,t]) 



(41) 
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It is readily observed that this operator depends on the curve s([0, 1]) C S in the parameter 
space S. One can treat it as an operator of parallel displacement along the curve s. 26-28 
For instance, if s([0, 1]) is a loop in S, the operator fl4"T| ) is the geometric Berry factor, and 
it can be treated holonomy control operator. ' 

APPENDIX A 

In order to complete the proof of Proposition |2|, let us first apply the relation ( p3|) to 
the Hamiltonian vector fields d a and i?j (p2|) . We obtain 

n a b = 5i # x a n\ = o, (42) 

tffflj = 5{, ^Q a k = 0. (43) 

The first of the equalities (|43"D shows that the matrix Q 3 k is nondegenerate, and so is i?*. 
Then the second one results in Qf, = 0. 

Using the well-known Kiinneth formula for the de Rham cohomology of a product, one 
can easily justify that the closed form Qz on W ([19]) is exact. Moreover, Qz = dE where 
3 takes the form 

Z = E a (J x ,y x )dJ a + E l (J x )d<p\ 

Of course, 3 is determined up to an exact form. Using the fact that components of oE = Jl^ 
are independent of y x and obey the equalities Q42|) - (|43]) , we obtain the following. 

(i) = — <9E a + <9 a 3j = 0. It follows that <9E a is independent of (p, i.e., 3 a is affine in 
V? and, consequently, is independent of <p since <p is a cyclic coordinate. Hence, <9 a Hj = 0, 
i.e., Hj is a function only of coordinates Jj. 

(ii) = — 5j5 +9 fe Hj. Similarly, one shows that S fe is independent of 9? and = <9 fc Sj, 
i.e., d k Si is a nondegenerate matrix. 

(iii) tit = -d b E a = 5%. Hence, E a = D a (J x ) + s a . 

(iv) fi^ = —dbE\ i.e., 3* is affine in s a . 

In view of items (i) - (iv), the Liouville form 5 reads 

5 = x a dJ a + [D\ J x ) + Bi( J x )s a ]dJi + Si( Jj) d^, 

where x a = s a + D a (J x ). Since the matrix <9 fc 3j is nondegenerate, one can introduce new 
coordinates Ij = 3j(Jj), 7 a = J a . We obtain 

3 = x<W a + [D H (I X ) + B':(I x )s a ]dI t + W- 

Finally, put 

0* = ^ - [i?«(J A ) + B^(J A )a«] 
in order to obtain the desired coordinates (p4|). 
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APPENDIX B 



In order to prove Proposition we first show that functions i^Fk make up a CIS on the 
symplectic leaf {V *Q, Qq) and iVo is its invariant manifold without critical points. Clearly, 
the functions igF k are in involution, and N is their connected invariant manifold. Let us 
show that the set of critical points of {i^F^} is nowhere dense in V *Q and Nq has none 
of these points. Let V *Q be equipped with some coordinates (q k ,p k ). Then the trivial 
bundle £ fll8|) is provided with the bundle coordinates (t,q k ,p k ) which play a role of the 
initial date coordinates on the momentum phase space V*Q. Written with respect to these 
coordinates, the first integrals F k become time-independent. It follows that 

dF k (y) = di* F k (£(y)) (44) 

for any point y G V*Q. In particular, if y G V *Q is a critical point of {igFfc}, then the 
trajectory £ -1 (?/o) is a critical set for the first integrals {F k }. The desired statement at once 
follows from this result. 

Since N is compact and regular, there is an open neighbourhood of N in V *Q iso- 
morphic to N x V where V C W 71 is a domain, and N x {v}, v G V, are also invariant 
manifolds in V *Q. 6 Then 

W" = C\N x V) = N x V (45) 

is an open neighbourhood in V*Q of the invariant manifold N foliated by invariant manifolds 
£,~ 1 {Nq x {v}), v e V, of the time-dependent CIS on V*Q. By virtue of the equality (||), 
the first integrals {F k } have no critical points in W" . For any real number r G (— e, e), let 
us consider a section 

h r : V*Q — ► T*Q, p o h r = -Hit, q j ,Pj) + r, 

of the affine bundle ( ([10|). Then the images h r (W") of H^" (|45| ) make up an open neigh- 
bourhood C/ of /i(AT) in T*Q. Because ((U) = W", the pull-backs (*F k of first integrals 
F k are free from critical points in U, and so is the function 7i* flip]). Since the coordinate 
r = p — h provides a trivialization of the affine bundle (, the open neighbourhood U of 
/i(iV) is diffeomorphic to the product 

h(W") x (-e,e) = /i(iV) x V x (-e,e) 

which is a trivialization of the fibration 

H* x (xC*F fc ) :E/-»Kx (-e,e). 

It remains to prove that the Hamiltonian vector fields of 7i* and (*F k on U are complete. 
It is readily observed that the Hamiltonian vector field 7t QITD of H* is tangent to the 
manifolds h r (W"), and is the image 7t = T/i r o j H o £ of the vector field 7// (fL4|). The 
latter is complete on W", and so is on U. Similarly, the Hamiltonian vector field 

Ik = -d t F k d° + d l F k di - diFkd' 
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of the function (*Fk on T*Q with respect to the Poisson bracket {,}t (0) is tangent to the 
manifolds h r {W"), and is the image 7/- = Th r o o ( of the Hamiltonian vector field 
of the first integral Ff. on W" with respect to the Poisson bracket {,}v flT2|). The vector 
fields on W" are vertical relative to the fibration W" — > R, and are tangent to compact 
manifolds. Therefore, they are complete, and so are the vector fields 7^ on U. Thus, U is 
the desired open neighbourhood of the invariant manifold h(N). 
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